home *** CD-ROM | disk | FTP | other *** search
/ The Complete Utilities To…ka 501 Killer Utilities! / 501 Killer Utilities! (Macworld July 1995).cdr / Programming / RLaB / examples / conv.r < prev    next >
Encoding:
Text File  |  1994-03-27  |  314 b   |  20 lines  |  [TEXT/RLAB]
  1. #
  2. # Simplistic convolution function.
  3. # x and y must be vectors...
  4. # No error checking.
  5. #
  6.  
  7. conv = function( x , y )
  8. {
  9.   local( X , Y , n, tmp );
  10.  
  11.   n = x.n + y.n - 1;
  12.   X = fft( x, n );
  13.   Y = fft( y, n );
  14.  
  15.   tmp = ifft( X .* Y );
  16.   if( type(x) == "real" && type(y) == "real") { tmp = real(tmp); }
  17.  
  18.   return tmp;
  19. };
  20.